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, In this work, we study the constraints on the scenario of light stop and sbottom in 
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£N| the next-to-minimal supersymmetric standard model (NMSSM), especially by a 125 

^"L GeV Higgs boson discovery and the LHC bounds on supersymmetry. The constraints 
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from dark matter detections are also taken into account. From the parameter scan, 
we find that the NMSSM can accommodate a light Higgs boson around 125 GeV and 



decay patterns well. The LHC direct SUSY searches with b-tagging are very powerful 

9-! 

JL and can set strong bounds on many NMSSM parameter points with light stop and 

(— 1 1 sbottom. We find i — > b\\ is a very promising channel for light stop detection if 

the mass splitting between xf an d Xi 1S very small. It is also pointed out that in 
order to close the parameter space of light stop and sbottom, new search strategies 
for signal channels such as pp — > ii^i — > tihhxiXi an d PP — > — > tiW + W~XiXi 



may be necessary. 

PACS numbers: 12.60.Jv,14.80.Ly 

I. INTRODUCTION 

The searches from ATLAS fl, CMS 2], and CDF as well as DO [g| established that 
there is a new particle around 125 ~ 127 GeV, of which its decay pattern is consistent with 
the predicted Higgs boson of the standard model (SM). If this new particle is a fundamental 
Higgs boson, it is reasonable to ask whether it is the predicted light CP-even Higgs boson 
in the supersymmetric models. However, a Higgs boson of 125 GeV seems a bit heavy for 
the minimal supersymmetric standard model (MSSM). In the MSSM, the tree- level mass of 
the lighter CP-even Higgs boson should always be smaller than the mass of Z boson, while 
the loop effects of stop can lift the Higgs boson mass up to 130 GeV or so. In order to avoid 
fine tuning problem, natural supersymmetry requires the third generation squarks are light 



Compared with the MSSM, a Higgs boson at 125 GeV can be realized more naturally 
in the NMSSM without confronting with the severe fine tuning problem, since it can have 
a much larger tree-level Higgs mass. This might be one of the reasons why NMSSM is 
appealing except that it solves the notorious \i problem in the MSSM. In Ref. 7j, the authors 
studied the constrained NMSSM with all parameters defined at the grand unification scale. 
It is found that the tension of the constrained MSSM for a 125 GeV Higgs boson mass and 
a light SUSY mass spectrum can be relaxed in NMSSM with more general parameters. 

Another reason that may favor NMSSM is the Higgs decay mode. The experimental 
data of the Higgs boson decay modes has shown a possible excess in the diphoton channel. 
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It may be a hint of new physics [8H10I]. In Ref. [Hi, H2J, it is shown that the NMSSM 
can accommodate the 125 GeV Higgs and enhance the branching of fraction BR(h — > 77) 
by reducing the branching fraction of BR(h — > bb) via mixing effect. A comprehensive 
study considering light third generation sparticles in the NMSSM also found the diphoton 
decay branching fraction can be enhanced [3]. Moreover, light charged sparticles such as 
charginos and charged Higgs boson may also contribute to the Higgs diphoton decay channel 
significantly 1^|. More studies on diphoton enhancement in both MSSM and NMSSM 
context can be found in 

LHC is searching for SUSY particles extensively in various channels. Up to now the null- 
result puts tension to many SUSY models. In Ref. 2l| it is found that the LHC sparticle 
search with 1 fb _1 at ^/s = 7 TeV exacerbated the tension between the 6(g — 2)^ USY anomaly 
and the branching fraction BR(B — > X s ^) in the constrained minimal supersymmetric stan- 
dard model (CMSSM). Of course, such conclusion is based on the assumed GUT relations 



among the soft breaking terms. When such relations are released as demonstrated in [22], 
the conclusion can be relaxed. For the pMSSM with 19-dimensional parameter space the 
current LHC search can not put very restrict constraints on supersymmetry in general. 

The LHC direct SUSY searches have set strong constraints on the masses of gluino and 
first two generations of squarks in mSUGAR and some simplified models 23|, |24] . Recently, 
there are many works on detection of light third generation squarks and gluino below 1 
TeV. It is pointed out that if the stop is light its si gna ture is possibly hidden from the 



LHC searches and can avoid severe 
scenario is still alive. More studies 
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27|. Therefore the natural SUSY 
28M35| are trying to improve the sensitivity of light stop 
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searches. For example, in 29|, |30|, the authors considered the stop dileptonic final state and 



explored a few kinematic observables to distinguish signal and background. The authors 



of Ref. 31J studied the top tagging technique for the stop search with semi-leptonic and 



dileptonic modes at the LHC. The hadronic top tagging technique has been examined in 



Ref. 32]. Interesting multiple lepton and jet final states from the multiple top decays are 
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. Especially 



investigated in [36[. Moreover, light sbottom searches are investigated in 
the sbottom-neutralino coannihilation scenario has been considered, it is found LHC has a 
good sensitivity to the parameter space by using tagged b-jet even for small mass splitting 
between sbottom and neutralino. 

Since in the NMSSM light stop and sbottom can be natural for a 125 GeV Higgs boson 
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131 ] . we did a systematic study to investigate the constraints on such a light stop/sbottom 
scenario from the LHC SUSY searches in this work. The constraints from B physics mea- 
surements and dark matter detections as well as the Higgs boson mass on the NMSSM 
parameter space are first considered. Then we study the constraints from various SUSY 
search channels at the LHC with ^ = 7 TeV and 2 ~ 5 fb _1 of data, including the jets + 
MET channel, associated monojet channel, lepton + jets + MET channel with and without 
tagged b jet(s). Since both stop and sbottom can decay into b jets, it is expected that b 
tagging should play an important role to distinguish signal and background. Our results con- 
firm this point and find the direct SUSY searches are powerful to exclude many parameter 
points with light stop and sbottom up to 500 GeV. 

The paper is organized as follows. In Section II, we briefly describe the NMSSM and our 
parameter scanning strategy. We will concentrate on the contribution of light stop to the 
mass of the discovered 125 GeV Higgs boson and the constraints from dark matter searches 
on the neutralino sector. In Section III, we analyze the LHC bounds from both ATLAS and 
CMS collaborations on the signatures of the light stop pair and sbottom pair production. 
Section IV is the discussions and conclusions. 

II. THE NMSSM 



In the NMSSM a singlet superfield 5* is introduced to solve the so-called ' 



The superpotential of NMSSM related to this singlet superfield S is given by 401 ] : 
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problem" . 



W NMSSM = XSH u Hd + + ■■■■> 
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where the dots denote the MSSM superpotential without /x term. When the electroweak 
symmetry is broken, the effective /i term can be naturally generated via the vacuum expec- 
tation value of S filed (labeled as vs), and can be written as /i e // = Xvs, which is expected 
to be of (9(100) GeV (of the same size as the rest of soft breaking terms). This /i e // can be 
traded off with m|, the mass parameter of the superfield S. The soft-breaking terms in the 
Higgs sector |40| is extended as 

Vnmssm = rh 2 H jH u \ 2 + rh 2 Hd \H d \ 2 + m 2 s \S\ 2 + (A x XSH u H d + ^A k kS 3 ) + h.c. (2) 

Compared with the MSSM, the Higgs sector becomes richer and contains three CP-even 
Higgs bosons, i.e. Hi, H 2 , and H 3 , and two CP-odd Higgs bosons, i.e. Ai and A 2 . Five 
new parameters A, k, A\, A k and n e ff are added compared with the MSSM. 



The Parameter space 



We use NMSSMTools 41] to perform a scan over the parameter space of the NMSSM. 
To obtain more generic conclusions, we consider the parameters defined at the electro-weak 
scale in our scan without assuming the unification of the NMSSM parameters at the GUT 
scale. We vary them in the ranges defined below: 

10" 4 < k < 0.5, 1 < tan (3 < 60, 50 < fj, < 500GeV, 

|A A |<4TeV, |A K |< 500GeV, lOGeV < M x < ITeV, lOOGeV < M 2 < ITeV, 
lOOGeV < m Q3 , m Ua < 2TeV, \A Us \ < 3TeV, lOOGeV < m t < ITeV. (3) 

There are a few comments in order on the ranges of NMSSM parameters. 

• The large A is helpful to raise the SM-like Higgs mass at tree-level and to ameliorate the 
fine-tuning issue confronted by the MSSM. When A tends to be zero, the singlet S will 
decouple from other Higgs fields. Under this limit, the phenomenology of the NMSSM 
may still be different from the MSSM due to the light singlet and singlino. These 
particles could affect the features of dark matter (DM). The decay modes of heavy 
sparticles produced at the colliders may also change, and some new search strategies 
will be necessary. Therefore we adopt two scan strategies which allow A varies in the 
ranges of [10 _3 ,0.1] and [0.1,0.8] with logarithmic and flat distribution respectively. 
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For the gluino and the first two generations of squarks, if their decay products are 
energetic jets and large MET can be reconstructed, the recent LHC results can put 
stringent limits on their masses in mSURGA and phenomenological SUSY, e.g. M3 ~ 



2 > 1.4TeV 23[ . If the first two generations of squarks are very heavy, the SUSY 



flavor and CP problems can be solved 42[. For the gluino, the naturalness of Higgs 



mass requires its mass should not be much larger than ~ ITeV js]. The main decay 
mode of gluino may be g — > tt/bb. For simplicity, in this work, we focus on the pair 
production of the third generation of squarks and leave this case for future study. 
Therefore we fix the soft-breaking parameters M3 = m^ 12 = 1.5TeV. To reduce the 
number of free parameters, we also assume the mass parameters of the third generation 
of right-handed squarks are the same, i.e. m£> 3 = mu 3 . 

In our scan, we require the mass of SM-like Higgs is in the range of 125 ± 2GeV and 
the SM-like Higgs boson can be either Hi or if 2- Except this, several phenomenology 
and astrophysics experimental limits are also considered. For flavor constraints, we 
require BR(B S -> X s7 ), BR(B+ t + za l )^R( j B X s fj.+/j,-), AM d and AM, sat- 



isfying experimental constraints at 2a 



43M45j . The theoretical uncertainties in these 



observables are considered as implemented in NMSSMTools. For BR(I? S — > fi + fi ), 
the upper-limits have evolved much in these two years. Recently ATLAS, CMS and 



,H 
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Cb collaborations have updated it to 2.2 x 10 -8 7.7 x 10 -9 (47j and 4.5 x 10 



-9 



at 95% confidence level, respectively, which are only a few times above the SM 
predictions. Here we adopt constraint BR(I? S — > fi + p~) < 4.5 x 10~ 9 given by LHCb. 
For the muon anomalous magnetic moment a^, we require SUSY effects can explain 
the discrepancy between SM prediction and experimental result at 2a. The mass lim- 
its for Higgs and charged SUSY particles from LEP, Tevatron and early LHC data are 
also adopted by in the NMSSMTools package [dl]. 

In our analysis, the lightest neutralino is required to be LSP and a candidate of the 
DM. Considering that there may be several types of DM in our universe, then the LSP 
in the NMSSM is just one kind of the DM, We only require the thermal abundance of 
neutralino satisfies a 3a upper-limit Q x h 2 < 0.1288 with the corrected DM relic density 
tth 2 = 0.112 ± 0.0056 reported by the WMAP 



49|. For this purpose, we define a 



fraction variable £ = Q x h 2 /Qh 2 , and the neutralino density in halo is p x = £,Pdm- The 
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DM detection limits given by experimental collaborations are obtained by assuming a 
certain DM density Pdm ~ 0.3 GeV cm -3 . Thus for the parameter points predicting 
£\/i 2 -C Qh 2 , these limits need to be rescaled. By doing this, we can examine whether 
the LSP in the NMSSM can sufficiently accommodate all data. The dark matter 
observations are calculated by MicrOmega 50] implemented in the NMSSMTools. 

Because this numerical scan is performed over a multi-dimensional parameter space, we 
use a Markov Chain Monte-Carlo (MCMC) method to increase the scan efficiency in our 
analysis. The total likelihood function L to t = n^Lj is evaluated by the likelihood functions 
based on the phenomenology and astrophysics experimental observables described above. 

i \2 



We define Lj = e 2<T * for two-sided constraints and Lj — 1/(1 + e CT * ) for upper-limits 



5l| , where Xi is the observable predicted by the model, /ijidj is the central values and error 



bars of experimental observables, a' taken as a' = 0.02/ij is the tolerance for upper-limit. 



B. 125GeV Higgs Boson 

In the MSSM, the tree level SM-like Higgs mass is smaller than Mz which is below the 
LEP limit m^ < 114 GeV. However, it can be lifted by the loop corrections (say top-stop 
corrections due the large Yukawa couplings). One loop formula for rrih is given by 

2 _ 2 i A 2 

m h — m h,tree + ^ m h,loop 

= Ml cos* 2 P + MLfMV4fl-4)l (4) 



t v t 



where v = 174 is the vacuum expectation value of SM Higgs, Mj = ^jm^m^ is related to the 
stop masses, X t = A t — fi cot is the stop mixing parameter. The Higgs mass depends on 
stop masses logarithmically. It means the stop bosons should be heavy in order to produce 
a large enough correction to the mass of the Higgs boson. We can also see the Higgs mass 
depends on stop mixing and is maximal for X = X t / Mi = \^6. 

In the NMSSM, the superpotential \SH u Hd can induce a new term X 2 v 2 sin 2/3 H u H<i in 
the Higgs potential. After rotating the 2x2 mass matrix of two CP-even neutral Higgs of 
doublets Hjj and Ho by angle (3, one diagonal element becomes Mf cos 2 2(3 + X 2 v 2 sin 2 2(3, 
which means that compared with the MSSM, the SM-like Higgs mass (mainly Hjj type Higgs 



boson) obtains a new contribution ~ A sin 2(3 



40|: 



ml = M 2 cos 2 2/3 + A t) sin 2 2(3 + Am 2 ^ + (5) 
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where dots denote the effects from mixing between Higgs doublets and the singlet. For large 
A with Xv > Mz, the is maximized for tan/3 = 1. Note that even sin 2/3 = 1 and A is 
taken as its maximal value ~ 0.7 which is consistent with the perturbation condition at the 
GUT scale, we still need a moderate stop loop contribution to obtain a 125 GeV Higgs. 




1000 
900 
800 




400 ff' 
300 
200 
100 

0.1 



0.4 0.5 

X 



w 



40 

30 ; 

20 
10 



FIG. 1: Left: mi versus X-t/Mi . Right: mj versus A. The color scale indicates A (left) and tan/3 
(right) respectively. 

It is interesting to ask how light the stop can reach after taking into account the 125 GeV 
Higgs boson in the NMSSM. To address such a question, below we show some scattering 
plots. In these plots, all the points have passed the constraints as being described in Sec. 

nm 

In the left panel of Fig. (TTJ, we show the correlation in the — X plane, the color 
bar shows the value of A. It is observed that for the large A ~ 0.6 — 0.7, the stop mixing 
parameter \X\ is allowed to be 1. When the light stop mass decreases, the stop mixing 
parameter will increase. If mi is below 300 GeV, X should be large than 3. It means there 
is a large mass splitting between two stop states. Because A can raise the mass of Higgs 
boson efficiently, it is easier to obtain a lighter t\ with a larger A. 

In the right panel of Fig. ([1]), we show the correlations between m~ tl and A, the color 
scale indicates tan /3. We also see that there are many points with small in the large A 
regime. As shown in Eq. [51 the Higgs mass depends on A sin 2/3, the values of tan (3 decrease 
at large A. We can see tan/3 < 30 for A > 0.3 and tan/3 < 10 for A > 0.6. 

Another interesting question is whether the NMSSM can accommodate the decay modes 



s 



and the decay branching fractions measured by the LHC collaborations appropriately. Espe- 
cially whether the NMSSM can explain the diphoton excesses and diboson (mainly h — > ZZ) 
data observed by LHC collaborations? Below we address this issue. 

In the SM, if the Higgs mass is determined, all the Higgs interaction couplings to SM 
particles can be obtained. In the new physics model, the Higgs couplings may differ to the 
SM predictions due to new parameters and particles. Therefore, the searches of Higgs partial 
widths at the LHC are very important to test the SM, and can provide crucial evidences of 
new physics. 

The effective Higgs couplings can be extracted from experimental data and can be com- 
pared with theoretical predictions. The parameter described the modifications of the Higgs 
couplings is defined as 

^hxx = ^hxxl^hxxi y°) 
where X can denote either heavy fermions, W boson, Z boson, photon or gluon. In the new 
physics model, both the production cross section and decay width of Higgs are rescaled by 
C\. The relevant Higgs partial widths would be determined by the ratio 

= ajpp -» h) NP BR(h -» XX) NP 
a{pp h) SM BR(h XX) SM 
For the 77, W + W~ and ZZ channel, the recorded Higgs events are mainly from gluon fusion 

gg —> h process. If the decay channel of 125 GeV Higgs is dominated by h — > bb, the Rhxx 
(X = W,Z,i) is approximated to be R hX x ~ C\ gg BR{h ->■ XX) NP /BR(h -> XX) SM ~ 
ClggClxxl^hbl- ^ or ^ ne ^ channel, additional electrons or muons are required to suppress 
huge QCD background, only the electroweak production channel qq — > hV is used to search 
Higgs signal. Therefore, the Rhxx can be given by R hb i ~ C\ vv BR{h — > bb)Np/BR(h — > 
bb) SM - 

In the NMSSM, the Higgs mass basis H™ ass = {Hi, H 2 , H 3 ] and interaction basis if* n * = 
{Ha, H u , S} are related by if™- - 55 = SiaH™ 1 . The reduced Higgs couplings to fermions and 
heavy gauge bosons can be given by 

S'i S- 

Chbb — Ch T f = 3, Chtt — —. — 01 Chvv — Sn cos/3 + S i2 sin/3. (8) 

cos p sin p 

In the SM, effective Higgs coupling to gluon Ch gg is dominantly determined by triangle top 
loop. In the SUSY model, stop loop would also contribute to Chgg-, it can be written as 

QSUSY _|_ qSUSY 

^hgg - J^sm + W> V*) 

^hgg, t 
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where Chxx, a is the loop contribution from particle A to the effective Higgs coupling Chxx, 
and C~ t is defined as C~ t = C*™[/C™ t . 

For the C^ 77 , the main contributions arise from W loop and top loop, which are related 

by C^w/Ci^t 8.3/1.8 in the SM. In the SUSY model, light charged particles (lj, 

such as light chargino, light charged Higgs boson, stop, sbottom, and stau, would provide 
additional contributions. We can also get C/j 77 approximately 



(jSUSY _|_ qSUSY _|_ qSUSY 

Chyy = hll '%sM h "?^SM — ~ — ~ l-28C hV v — 0.28(CW + Q) + Cf + C%+, (10) 



where we have used the relations C^ S1 ^/C^ W ~ C hVV , C^ s J/C^ >t ~ C htl and 
^h US } ' /^hrrt t ~ ^hg 9 S J l^hgg.t (neglecting the high order QCD corrections induced by 
stop), and Cf is defined as Cf = C^ g s J /C^. Note that Cf is proportional to Yukawa 
coupling yh Tf and inversely proportional to stau mass m T , it means Cf is important for large 
tan (3 > 50 ~ 60 and light stau 8|. 

An important feature of Higgs phenomenology in the NMSSM is the exotic Higgs decay 



modes to light neutralinos or scalars 



52 



53J . If Hi is SM-like, H 2 might be much heavier than 



125 GeV. If Hi is SM-like, H\ and A\ would be light singlets. In this case, the branching 
ratios of H 2 — > XiXii H 2 —> H\H\ and H 2 — > AiA± could be large due to Higgs mixing 
and kinematics. For the invisible Higgs decay H 2 — > XiXi: the possible search channels 
are // — > Vh — > leptons+MET or gg — > hj — > mono-jet+MET. For the decay channels 
H 2 — > H\H\ and H 2 — > A±Ai, the light singlets would decay into rf or bb. The final 
states are four fermions which can be significant for some parameter points. However, these 
exotic decay modes might suppress the standard Higgs decay modes to heavy gauge bosons, 
photons and bb, and could be tested by the global-fitting of Higgs decay partial widths in 
the future. 

To examine the question whether the NMSSM can accommodate the diphoton excess, we 
show scattering plots to demonstrate the correlations between the effective couplings and 
i? 77 . In the left /middle/right panel of Fig. (T5]), we demonstrate the correlations between 
i? 77 with Chbi/Ch-yy/Chgg, the blue/red points denote SM-like Higgs is H\jH 2 . We can see 
CVy 7 and Chgg always vary in the range of ~ 0.8 — 1.1, and _R 77 is sensitive to C hb i. We also 
find i? 77 is inversely proportional to R b i and R h i is larger than 1 for i? 77 > 1. 

These results can be understood in term of explanation given above. In the decoupling 
limit Ma 3> Mz {Ma = is the mass of heavy CP-odd Higgs), the main component of SM-like 
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FIG. 2: The correlations between i? 77 and C hbb /Ch- ry /Ch gg are demonstrated. The blue/red points 
correspond to H1/H2 is the SM-like Higgs boson. 



Higgs is H u due to S 1 ^ ~ sin/3 ~ 1 for tan/3 ^> 1. From Eq. (jSJ), we can get Chit and Chww- 
If the stop mass parameter -Jm^m^ is much larger than m t , the stop loop would not provide 
larger contributions to Ch gg - Then Ch gg ~ 1 and C^ 77 ~ 1 can be obtained from Eq. (JUJ) and 
( TTU1) respectively. Depending on parameters in Higgs sector, can be significant reduced 
by mixing effect. As we have mentioned, if the main SM-like Higgs decay mode is h — > bb, 
and R hX x(X = W,Z,j) is approximately to be R h xx ~ c l gg c lxx/ c l h i ~ l / C h\v This 
relation explains the inversely proportional correlation between i? 77 and C hbb f° r > 0-8 
as shown in Fig. (j2j). This is the case for H\ is SM-like Higgs. 

If H2 is SM-like, the branching ratios of exotic decay modes might be significant, the i? 77 
would be much suppressed even the C h bi is still ~ 1- For either blue/red points, to decrease 
R h l and increase i? 77 can be done by tuning the singlet-doublet mixing parameters. 

We also show the correlations between i? 77 and Ryv i n Fig- ©> the color bar indicates 
C h bi {Rbt>) m the left (right) panel of Fig. ()3]). From the left panel, we can also see the R b i 
is smaller than 1 for i? 77 > 1. Moreover, the i? 77 is always proportional to Ryy due to the 
approximations Rhw/R-yy ~ ^hw/^h-ry an< ^ ^77 ~ l-28CVv — 0.28(7^+ Csusy- As we have 
mentioned above, R h bi depends on Cyv due to the production process. As Cvv is almost 
~ 0.8 ~ 1.1 shown in Fig. ([2]), R hb i is mainly determined by BR(/i — > bb). For the large 
i? 77 > 1, BR(h — > bb) is suppressed due to small C hb i. We also find the i? 77 can be smaller 
than 0.5 due to the significant Higgs exotic decays, both the Ryy and R bb are suppressed in 
this case too. 

According to the analysis given by CMS group, the ratio of the couplings of Higgs to 
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FIG. 3: The correlations between Ryv, and C b i/R b i are displayed in the left/right panel. The 
observed values given by ATLAS and CMS are also shown. 

fermions is around 0.5 ± 0.3. In the NMSSM, this suppression can be accommodated by 
the mixing between the singlet and doublet Higgs boson, while keep the couplings to vector 
weak bosons close to one. In Fig. ([3]), we also mark out two current global values of Rzz 
(0.7 ± 0.4), ~ (1.6 ± 0.4) and R zz ~ (1.4 ± 0.5), R yy ~ (1.8 ± 0.4) given by CMS 
and ATLAS |l| respectively. For the light stau region which may be helpful to ease this 
tension, we find the 5(g — 2)^ and flavor physics put stringent bounds to the stau mass 
and tan/3 in our scanning. The parameter region providing the enhancement by the light 
chargino/charged Higgs boson to i? 77 has not been reached since the parameter A is confined 
to be less than one in our scanning. 

It should be noticed that with the current experimental error bars and statistics, it is too 
early to conclude that the decay patterns of Higgs boson have confirmed the existence of 



new physics. As pointed out in Ref. [54j, the large uncertainty in the parton distribution 
function can also affect these results. Future data and reduction in the uncertainty are 
needed to make sure whether the new physics has been indicated in the Higgs decay modes 
already. 



12 



C. Dark Matter Bounds 

In this section, we will discuss the constraints from the DM detections. We assume the 
LSP and DM candidate is the lightest neutralino. Because the neutralinos have additional 
singlino component in the NMSSM, the phenomenology of DM is different from that in the 
MSSM. Especially, the LSP can be pure singlino. In this case, the LSP can be lighter than 
100 GeV and can easily escape the constraints from the invisible Z decay measurements due 
to its almost vanishing coupling to the Z boson. Below we will address the issue whether 
the singlino in the NMSSM can help to ease the tension between the theories and the 
experiments. 
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FIG. 4: The Singlino/Bino/(Wino+Higgsino) content of the LSP versus LSP mass. The color scale 
indicates the LSP relic density. 

First, Let's examine the constraints from the dark matter relics density. As we know, the 
neutralino mass basis X° =1 _ 5 and interaction basis ip® = {B, W, H d , H u , S} are related by 
X° = Nijipj. In the left /middle/right panel of Fig. (j3J), we show the singlino/bino/Higgsino 
and wino content |iVi5| 2 ( | iVi 2 1 2 + | ^V^i 3 1 2 + |A r i4| 2 ) of the lightest neutralino, the color scale 
indicates neutralino thermal relic density Q x h 2 . To obtain the suitable relic density Q x h 2 < 
0.1388, the neutralinos should have enough annihilation cross section (av) > 3 x 10 _26 cm 3 s _1 . 
In the left panel, we can see the most of lightest neutralinos with significant singlino contents 
are lighter than 100 GeV. Because the mixing terms between singlino and Higgsinos in the 
mass matrix are proportional to A, the LSP as pure singlino means A is small. In this case, 
the singlino mass is approximately 2nfi/\, the lightest CP-even Higgs and CP-odd Higgs 
are also singlet-like with small masses. The main DM production mechanism might be 
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annihilation with s-channel Z resonance or CP-odd/CP-even Higgs resonance 1 . This is the 
reason why there are many points condensing around the range of ~ 40 — 70 GeV. Moreover, 
if the A is not very small, neutralinos could annihilate into light scalar pairs H\H\, A\A\ or 
H\A\ via the t-channel by the x° exchange or the sufficient large singlet- singlino interaction. 

If the lightest neutralino has sufficient Higgsino or wino content, the neutralino annihila- 
tion cross section can be large via t-channel chargino exchange to W + W~. This is the case 
for points with Q x h 2 <C 0.1 and m x > 80 GeV in the right panel of Fig. (j3J). Because the 
lightest chargino xt can be pure wino or Higgsino. It is also possible to find the parameter 
points with almost degenerate neutralino Xi an d chargino xt which means neutralino can 
obtain suitable relic density via large co-annihilation XiXi an d XiXi ■ 

Next we consider the constraints from the indirect astrophysics search. If the DM particles 
annihilate into heavy quarks, charged fermions and gauge bosons at the present time, these 
annihilation final states can induce significant gamma ray flux which can be detected by 
air shower cherenkov detectors or satellite detectors. In Fig. ([5]), we plot parameter points 
in the m x o vs. £, 2 (crv) plane. The color scale indicates Q x h 2 . Note that the DM thermal 
averaged annihilation cross section is rescaled by £ 2 because the gamma ray flux depends 
on DM density square and the actual neutralino density may be ^Pdm- in Fig. fl5]), the 
upper-limits derived from Fermi gamma ray observations towards dwarf spheroidal galaxies 



are also shown [56]. These limits are combined by null results from ten dwarf galaxies, 
and only a few times above the "natural value" (av) = 3 x 10~ 26 cm 3 s~ 1 for DM with 
uidm ~ O(10 2 ) GeV. For parameter points with large wino and Higgsino contents and then 
large (av) > 3 x 10~ 26 cm 3 s _1 , the thermal neutralino density ^Pdm is very small. Thus 
the reduced annihilation cross section rescaled by £ 2 can escape the constraints easily. The 
gamma ray limits become more stringent when DM mass decreases. However, for light 
neutralino m x o < 70 GeV with significant singlino content, the annihilation process via s- 
channel Z or CP-even Higgs exchange is p-wave which is much suppressed at the present 
universe. Therefore we can see these limits do not exclude many parameter points 2 . The 



1 For the s-channel Z exchange annihilation, small Higgsino content is still needed because ZxiXi coupling 
is proportional to |7Vi 3 1 2 - |iV 14 1 2 . 

2 The (av) is often estimated at v ~ 1CP 3 which is the typical DM velocity in the present Galactic halo. 
However, the velocities of DM particles in the dwarf galaxies are about an order of magnitude smaller. 
When the gamma ray limits from dwarf galaxies are taken into account, this effect needs to be considered 
for the velocity dependent annihilation cross section (55| . This effect may enhance the cross section signif- 
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more stringent limits can be derived by Fermi gamma ray observations towards Galactic 
center, because DM particles are more condensate in this regime. The main problem is how 
to subtract the complicated astrophysical backgrounds precisely. However, it is possible to 
improve the constraints to be below the "nature value" for O(10 2 ) GeV DM in the future. 




FIG. 5: ^ 2 av versus m^o. The color scale indicates the neutralino relic density. The black triangles 
are the LSP with singlino component dominated ( | TVi 5 1 2 > 0.8). The upper limits on three annihi- 
lation channels XX ~ > bb, tt and W + W~ given by Fermi-Lat dwarf galaxies observations are also 
shown 56]. 



Then we examine the constraints from the direct searches. We consider both spin- 
independent and spin- dependent cases. For the spin-independent constraints, we focus on 
the constraint from the XENON100, which is the most stringent bound for dark matter di- 
rect searches. For the spin-dependent constraints, we include bounds from the neutrino flux 
measurements. In Fig. fl6]), we plot parameter points in the m x o vs. £,<Jsi plane. The color 
scale indicates Q x h 2 . Here £,<Jsi is the reduced spin-independent neutralino-nucleon elastic 
scattering cross section. We also show the most stringent constraints set by XENON 100 



in 2011 



57] and 2012 581 ] . We can see the XENON limits have excluded many parameter 



points with corrected DM relic density Q x h 2 ~ 0.11. The spin-independent scattering pro- 



icantly when the main annihilation process is s-channel CP-odd Higgs exchange. The large enhancement 
often requires a very narrow Breit-Wigner resonance with tiny mass splitting parameter 1 1 — m\/ 4m^ | <C 1 
and decay width F^/m^ <C 1. This effect would not change our results very much, and is neglected here. 
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cesses via squark exchanges are strongly suppressed due to the heavy squarks assumed to 
avoid collider constraints. The main process is the exchange of Higgs through t-channel. 
The cross section of such process depends on the wino and Higgsino contents of the lightest 
neutralino and the Higgs masses. Therefore, the neutralinos with masses of m x > 80 GeV 
and intermediate wino and Higgsino contents are strongly disfavored by direct detections, 
unless the agi is reduced by small £. Moreover, if the neutralinos have significant singlino 
contents, and lightest CP-even Higgs is singlet-like, the neutralino-nucleon scattering can be 
enhanced by additional Higgs-neutralino couplings and small Higgs mass in the propagator. 
These parameter points might also be easily excluded by XENON limits. 




FIG. 6: £<7s7 versus m^o. The color scale indicates the neutralino relic density. The black tri- 



angles are the LSP wit 



XENON 100 in 2011 



i singlino component dominated ( | iVi 5 1 2 > 0.8). The upper limits set by 



57J and 2012 



581 ] are also shown. 



For the spin-dependent neutralino-nucleon scattering, the constraints established by di 



rect detections are very weak. The most strong direct constraints given by Coupp 
KIMS 



and 



60| are of the order of O(10 x ) pb. More stringent limits can be set by high energy 



neutrino telescopes. If the DM particles lose their energies by scattering with solar nucleons 
and are trapped in the center of the sun, they could annihilate into to heavy fermions or 
gauge bosons and produce detectable high energy neutrino signatures. The signature flux 
can be determined by DM-nucleon spin-dependent scattering rate. Moreover, it also de- 
pends on the fractions of certain DM annihilation channels /, = cr(xx ~^ XiXi) I a xx w hich 
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could produce neutrinos. In the left (right) panel of Fig. [7J we present the parameter points 
in the m^o vs. fen (m^o vs. (,o- SD f vv ) plane. The limits given by Coupp |59|, Superk 
6l| and IceCube 62] are also shown. Note that the neutrino energy spectra from different 
DM annihilation channels are different, the limits from neutrino telescopes are derived from 
experimental results for different channels. Because the neutrino energy spectra from bb (qq, 
rf) are much softer than those from W + W~ (ZZ, tt), the limits from W + W~ channel are 
more stringent. From the left panel of Fig. ([7]), we can see recent direct detections and 
neutrino signatures from soft channels do not constrain models very much. However, if the 
neutralinos have intermediate Higgsino contents, the spin-dependent cross section might be 
large due to ZxiXi coupling which is proportional to | iV"i3 1 2 — | A^i4 1 2 . In this case, some 
parameter points with £,<Jsd > O(10~ 4 )pb and fyy ~ 1 have been excluded. The expected 
limit which can be set by IceCube 86 strings are also shown. We can see the IceCube has the 
capability to exclude most of the parameter points which predict neutralinos with significant 
Higgsino contents and correct DM relic density Q x h 2 ~ 0.11. 
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FIG. 7: Left: £,<jsd versus m^o. Right: £o-sDf(XiXi ~^ VV) versus ~m^o. The color scale indicates 
the neutralino relic density. The black triangles are the LSP with singlino component dominated 



(|iV 15 1 2 > 0.8).. The limits given by Coupp 



59(], Superk 



6ll ] and IceCube 



62 ] are also shown. 



To examine the impact of singlino to the dark matter searches, in Fig. (IMZI), we have 
marked out those allowed points with dominant singlino constituent by the black triangles. 
It is observed that when the LSP is singlino dominant, it can escape the constraint from 
Fermion LAT easily, as demonstrated in Fig. ([5]). While the bounds from XENON100 2011 
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and 2012 are really impressive and can exclude a certain fraction of those allowed points 
even the LSP is singlino dominant. For the neutrino bounds, since many LSP with singlino 
dominant are lighter than W boson, they can be still consistent with data due to weak 
constraints from bb channel. For the LSP heavier than W boson, the annihilation cross 
section of x°X° ~^ W + W~ is always small enough and can be safe. From the analysis shown 
above, we can see that the singlino/bino in the NMSSM can help to ease the tension between 
theories and experiments, while the wino and Higgsino like dark matter candidates are more 
constrained. 



III. LHC SUSY SEARCH BOUNDS 
A. The production and decay of light stop/sbottom at LHC 

In this section we study the SUSY search bounds on the light stop/sbottom pair signatures 
from LHC. For our purpose below we make two additional requirements for the parameter 
points past all constraints in Sec II. First, we require the stop mass should be lighter than 500 
GeV in order to have a large enough production rate. Since we have assumed mu 3 = m£> 3 , so 
the lightest sbottom is light in our analysis. Then the largest color sparticle signatures are 
the stop pair and sbottom pair productions. Second, to accommodate the SM-like Higgs in 
the NMSSM, we require Ryv + -R77 > 1-8 in the parameter space. Implicitly, this condition 
also means that the branching ratios of the exotic Higgs decay modes should not be very 
large. Consequently, it also suggests that the lightest neutralino, CP-odd and CP-even Higgs 
can not be very light. This feature will also affect the decays of other sparticles. After taking 
into account these two extra requirements, we choose 647 parameter points that survived all 
our criteria for our simulations. Below we study the constraints from direct SUSY searches 
at LHC to these points. 

In the left panel of Fig .([8]), we show the cross section of s top pair production an at the 



LHC with yfs = 7 TeV. Here we use the package Prospino2 [63| to calculate a% including 
the next leading order corrections. Because the main production channel of stop pair at the 
LHC is gg — > tt, we can see an is uniquely determined by stop mass. For the stop with mass 
of m t ~ < 500 GeV, the an is larger than 45 fb, and there would be more than ~ 200 stop 
pair events at the LHC with 5 fb -1 of data. 
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FIG. 8: Left: the NLO cross section of stop pair production as a function of stop mass. Right: 
sbottom mass versus stop mass, where the color scale indicates neutralino mass. 



The masses of stop and sbottom are shown in the right panel of Fig. (jSJ) where the 
color scale indicates the lightest neutralino mass. We can see the mass splitting between 
stops and sbottoms are small due to our assumption mu z = rriD 3 and the lighter stop and 
sbottom quarks are either dominantly left-handed or right-handed. It is supposed that this 



parameter configuration can pass the electroweak precision tests easily [13|, \64 \. 




100 150 200 250 300 350 400 450 500 100 150 200 250 300 350 400 450 500 
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FIG. 9: The lightest neutralino/chargino mass versus stop mass in the left/right panel. The color 
scale indicates the branching ratio of stop BR(i — > txi) /BR(t — > bxf) in the left/right panel. 

The stop decay pattern is dominantly determined by the mass splitting between the stop 
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and light neutralinos/charginos. If the stop is much heavier than the lightest neutralino 
and chargino, the main decay channels are two body decays t — > tXi and t — > b\i ■ In Fig. 
([9]), we show the relations between the sparticle mass spectra and the branching ratios of 
t — > t\\ and t — > bxf . We can see the decay modes depend on the neutralino and chargino 
mass spectra. The processes of stop decay into txlL ^xi an d bxt might be also significant for 
heavy stop (as demonstrated in first two bench mark points given in Table. (1VIII) ). In this 
case, even the LSP is very light < lOOGeV, decay mode t — )■ t\\ may also be suppressed. 
These processes have longer decay chains than t — > txi, and produce softer final states. 
When the mass splitting is too small to forbid above two body decays, the three body decay 
channels t — > bW\\ and t — > bvljblv become important. If these processes are also kinematic 
forbidden, the loop induced two body FCNC decay t — > cxi would be dominant. Moreover, 
the four final state decay modes t — > &J1J2X1 an d t — > Mu^Xi are also possible if the W boson 
in the three body decay mode is not on-shell. 

The main two-decay modes of sbottom include b\ — > bx® and b\ — > tx~ ■ The three-body 
decay modes include b\ — > t*x~ — > bW + x[ ■ The decay chain of sbottom can be quite long 
and all the final state from its decay can be soft as demonstrated in the third bench mark 
points in Table. fIVIIj) . 



B. constraints on stop/sbottom pair signatures 

The general SUSY search bounds for squarks and gluino have been provided by ATLAS 
and CMS, corresponding to an integrated luminosity of 2 ~ 5 fb -1 . In this subsection, 
we investigate the constraints on light stop pair and sbottom pair productions based on 
these results. In our studies, parton-level events pp — > tt and pp tt + jets are generated 



by MADGRAPH5 [65]. PYTHIA 66[ is used to perform the parton shower, decay, final 



state radiation, and hadronization processes. The detector effects are simulated by PGS4 



67|. To avoid the double-counting issue, we adopt the MLM matching scheme and choose 
Qcut — 80 GeV in our simulation. Jet candidates are reconstructed by using the anti-kt jet 
algorithm (which is infrared and collinear safe) with a distance parameter of 0.4/0.5 for the 



ATLAS/CMS searches 



3 



3 The basic selected conditions for jets and leptons are slightly changed in different searches. In general, 
these conditions are p t > 20-40 GeV, | ^7 1 < 2.5-3 for jets, p t > 10-25 GeV, \r]\ < 2.0-2.5 for electrons 
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Currently, most of experimental groups from both CMS and ATLAS collaborations work 
in the simplified model. The signals are assumed to be pp —> bib\ —> bbXiXi or PP ~~ t\t\ ~^ 
tt-XiXv It is interesting to examine what might happen in a concrete model like the NMSSM 
and how new decay modes can be affected by these direct searches. We will consider two 
categories of constraints from the direct search: 1) the searches for the final states without 
b-jets; 2) the searches for the final states with b-jets. 

1. Constraints for Final states without b-jets 

If the dominated decay modes of squarks and gluino are q — > qxl and g —> qqxl, the 
main features of events are energetic jets and large MET. This signature channel can set 
the most stringent constraints on CMSSM and simplified model without light stop/sbottom. 
For the light squarks of the third generation, such signatures would be suppressed due to 
smaller production cross section and different decay modes. Because such analysis requires 
very hard jets and large MET, the events of stop pair with many soft jets can escape the 
constraints easily. Below we investigate how the m e ff and -^cuts as well as the associated 
mono-jet searches can affect our selected points. 



and muons. Moreover, the electron candidates in the barrel-endcap transition region with 1.44 < rj < 1.57 
are rejected. Here we used the basic selections as adopted by ATLAS and CMS according to different 
researches, and do not list them in the following discussions. 
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Requirements 


A 


A' 


B 


C 


D 


E 


pr [GeV] > 


160 


N Jet (p T > 130 GeV) > 


1 


Nj e t{pT > 60 GeV) > 


2 


2 


3 


4 


4 


4 


NjetiPT > 40 GeV) > 










5 


6 


&<l>{jii$r)rnin > 


0.4 (i=l,2,(3)) 


0.4 (i=l,2,3), 0.2 {jP T > 40 GeV) 


pr/m eff {N Jet ) > 


0.3 (2j) 


0.4 (2j) 


0.25 (3j) 


0.25 (4j) 


0.2 (5j) 


0.15 (6j) 


m e //(incl.) [GeV] > 


1900/1400/- 


-/1200/- 


1900/-/- 


1500/1200/900 


1500/-/- 


1400/1200/900 


N obs 
hm 


2.9/25/- 


-/29/- 


3.1/-/- 


16/18/58 


10/-/- 


12/12/84 



TABLE I: Summary of cuts and observed 95% CL upper limits on the excess event number, 
following the ATLAS jets+MET analysis for 4.7 fb -1 ^]. 



The first constraint is from the jets plus missing energy searches. We list the cut condi- 
tions adopted by ATLAS collaboration in Table. ()T]) 23J- This analysis is based on 4.7 fb _1 
of data, and all the events with isolated electrons or muons are rejected. The azimuthal 
angle A0(jj,^) is defined as the azimuthal angle separation between the fir and the jets. 
The effective mass is defined as 

Nj N, 

m eff = fir + J2PT + J^PT- (11) 

i=l i=l 

It is obvious to see the effective mass characterizes the mass scale of SUSY particles directly 
produced by pp collisions. Large m e ff ~ ITeV is helpful to reduce the SM backgrounds, 
such as W + jets, Z + jets, ti and single top, but it also suppresses light stop/sbottom pair 
events with soft jets severely. 

For each signal region, three m e // cut conditions (denoted by 'tight/medium/loose') are 
taken into account. In Table. ([!]), the 95% CL. observed upper limits Nu m on number of 



new physics events given by experimental collaboration are also listed |23|. For comparison, 
we show the ratio of predicted event number to upper limit Nu m in Fig. (jTt)]) where 
the color scale denotes the LSP mass 4 . Here we have summed the signature of stop pair 
and sbottom pair together. If the m e // cut is chosen to be too large, most of the signals 



4 It should be noticed that the bounds shown in our analysis have not included theoretical and experimental 
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ATLAS jets+MET 4.7fb" Region A', medium 







250 


10 f 




- 


200 








150 


e 1 : 


= 




100 


Z " 0.1 r 






50 


0.01 - 



100 150 200 250 300 350 400 450 500 

m,(GeV) 



i 



100 150 200 250 300 350 400 450 500 
m,(GeV) 



ATLAS jets+MET 4.7fb" 1 Region C, medium 



.V. - >itf.v« :>.*.: 
.a .•.*-'%t'*SV* 



n 



i 



ATLAS jets+MET 4.7fb Region E, medium 



250 
200 



100 z 



100 150 200 250 300 350 400 450 500 
mj (GeV) 




100 150 200 250 300 350 400 450 500 
m,(GeV) 



FIG. 10: The ratio of tt*+bb* event number and observed 95% CL upper limit versus stop mass in 
the signal regions A, A', D, E. The color scale indicates neutralino mass. 



would be rejected. Otherwise, there are too many background events which leads to a much 
weaker the upper-limits for new physics. Better constraints can be demonstrated by the 
four medium channels A, A', C and E in Fig. ffTO]) . 

We have also checked the results for all the channels, and find jets+MET channel can set 
marginal constraints on the most of parameter points when the stop mass is below ~ 200 
GeV and the main decay mode is t — > c%°- The large cross section of stop pair production ~ 
O(10) pb can induce event excess, as indicated in A and A' region in Fig. (IIU|) . Interestingly, 
it is noticed that some parameter points with light LSP and stop can not be excluded by this 
search. These parameter points also predict light charginos, sleptons or second neutralino, 
therefore the decay chain can be long and the final states contain less hard jets, which can 
hide into the background events, as shown in all regions of Fig. (I10p . 



errors. Although our simulated results are close to those experiment ones, it should be remembered that 
when these uncertainties are taken into account, the lines should become to be bands. 
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We can see there are quite a fraction of points in A and A' situated below N/Nu m = 0.1, 
while all those points are above N/Nu m = 0.1 in C and E. For the heavier stop with 
dominated decay mode t —> txi, the channels requiring high jet multiplicity should be 
more efficient. For instance, we can read out the values of N/Nu m is only close to 1.0 for 
red parameter points with m 4 - > 350 GeV and m^o < 100 GeV in the region E-medium 
presented in Fig. (flO]) . 

We also check the constraints from ATLAS jets+MET research based on 1 fb _1 of data 
. In this analysis, the m e ff cut is required to be O(100) GeV. Considering the huge SM 



backgrounds prevent from setting a better constraint to stop/sbottom mass, we find the 
N/Num is much lower, e.g. <O(10~ 1 ). 



N Jet (p T > 100) 


Nj e t{pT > 30) 




fir 




> 1 


< 2 


< 2.5 


> 250/ 300/ 350/ 400 


600/ 368/ 158 /95 



TABLE II: Summary of cuts and observed 



CL upper limits on the excess event number, 



following the CMS mono-jet+MET analysis for 4.98 fb" 1 



The second constraint which will be considered here is from the associated mono-jet 
production. As is well-known, when the dark matter particles are directly produced by pp 
collisions, one possible search channel is the mono-jet + MET 69]. The mono-jet is produced 
by the initial state radiation and can be energetic. If stop is nearly degenerate with the LSP, 
the soft jets from stop decay might not be reconstructed by the detectors. In this case, the 



stop production ttj can be constrained by mono-jet + MET research [27. l28l . I70l|. The cut 
conditions and upper-limits given by CMS are summarized in Table. (JXXJ) 71]. This analysis 
is based on 4.98 fb" 1 of data, and all the events with isolated electrons or muons are rejected. 

In Fig. (HID, we show the ratio N/N Km with both fi> 350 GeV and fi> 400 GeV cases. 
For these large fir cut conditions, the jet from the initial state radiation is required to be 
very energetic pj ~ fir- It is obvious to see only light stop with rrif < 200 GeV and large 



production cross section can have large event rate, as demonstrated in Fig. (llll) . For heavier 
stops with t — > tx^l, cut condition on the third jet will suppress the events with high jet 
multiplicities and leads to a weaker constraint. For the sbottom which is not degenerate 
with LSP, the dominated decay mode can be b — > bxi- Therefore, the cuts on transverse 
momentum of jets can be satisfied easily. But the sbottom pair can not induce large fir 
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signal as required in this analysis which also leads to a weaker constraint. 
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FIG. 11: The ratio of tt*+bb* event number and observed 95% CL upper limit versus stop mass 
in the signal regions for f}p > 350 GeV and ftp > 400 GeV. The color scale indicates neutralino 



mass. 



2. Constraints for Final states with b-jets 

In this subsection, we explore the impact of b-jets+MET searches on stop/sbottom pair 
production signatures. The identification of b-jet is helpful to reduce the huge QCD back- 
grounds. In this case, the dominated SM backgrounds are top pair production and associated 
production of W/Z with heavy flavor jets. The di-bosons production WW, ZZ and WZ 
are sub-dominated due to smaller electro-weak cross section. 



As discussed in Ref. jfj|, the "naturalness" of SM-like Higgs mass suggests a light stop 
< 700 GeV and a not very heavy gluino ~ 1 TeV in the SUSY model. In this case, 
the gluino pair production has a moderate cross section, and the cascade decay productions 
would contain many top and bottom quarks. For example, the typical SUSY search channels 
are gg -> tttt* -> tttt + MET or gg -> bbbb* ->■ bbbb + MET. The multi-b jets in the final 
states are very powerful to suppress the SM backgrounds. Therefore, the LHC has strong 
capability to test or exclude such scenario. 

In this work, we have assume the gluino is very heavy > 1.5 TeV. Therefore the main 
production signatures are pp — > tit* and pp — > b\b\. For the light sbottom pair production, 
the b-jets+MET search can constrain the channel bb — bb + MET if the Am~ b = m~ b — m^o 
is large enough (say larger than 50 GeV). For the stop pair production, if the dominated 
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stop decay mode is t — > bWx or t — > tx — > bWx, the b-jets in final states are less energetic. 
To pass the p T cut on the leading b-jet, the Am ( - = m 4 ~ — m^o is required to be large, and 
the detectable capability is lower than that of sbottom pair. Here we point out that if the 
chargino is light and light stop has large left-handed component, the t — > b%i can become 
significant. Especially, if the chargino mass is nearly degenerate with LSP which can be used 
to obtain required DM relic density through co-annihilation mechanism, the kinematics of 
t — > bxt is very similar to the b — > bxi- In this case, the b-jets+MET search is also useful to 
test or exclude stop pair signatures. This feature is also addressed in the first bench mark 
point in Ref. fl- 



bjeti 
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j4 3 
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TABLE III: Summary of cuts and observed 95% CL upper limits on the excess event number, 
following the ATLAS 2b-jets+MET analysis for 2.05 fb _1 [73|. 
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FIG. 12: The ratio of tt*+bb* event number and observed 95% CL upper limit versus stop mass in 
the signal regions for mcT > 100 GeV and mcT > 200 GeV. The color scale indicates neutralino 
mass. 

First, we consider the ATLAS 2b-jets+MET search based on 2.05 fb _1 of data [3]. This 
search is optimized for sbottom pair production with sbottom branching ratio of BR(b — > 
bxi) = 100%. The corresponding cut conditions and limits are summarized in Table. (1IIID . 
In this analysis, the number of jets with p t > 50 GeV is required to be exactly two. No m e ff 
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ATLAS 2bjets+MET 2.05fb'\ m CT >100GeV ATLAS 2bjets+MET 2.05tb"' , m CT >200GeV 




FIG. 13: The ratio of tt* event number and observed 95% CL upper limit versus stop mass in the 
signal regions for mcT > 100 GeV and mcT > 200 GeV. The color scale indicates chargino mass. 

cut is imposed because sbottom should be light to provide large production cross section. 
A boost-corrected con-transverse mass ttict is introduced \l\ 

mcT = v/(^+^) 2 -(^-^) 2 , (12) 

where ttict is invariant under contra-linear equal magnitude boosts. It can be expected that 
the distribution of rncr display an endpoint at (m~ — m? )/^ when two b-jets are co-linear. 
In Fig. (I12jl . the ratios N/Nu m are shown . We find this search can exclude many selected 
parameter points, especially when the mass splitting Am^ = — m^o is large enough with 
rrit in the range from 250 GeV to 350 GeV. When comparing the left and right panel of 
Fig. (fT2|) . it is easy to read out the fact that the smaller mcr is more sensitive to smaller 
stop/sbottom mass region. It might be interesting to notice that the maximum excluded 
mjj can reach to 380 GeV by this observable. 

Moreover, we show the ratios N/Nu m only for stop pair production in Fig. ( {TBI where 
the color scale indicates the mass of chargino Xi ■ It is noticed that even if sbottom is very 
heavy, the 2-bjets search can be useful to put constraints on stop pair production. This can 
occur for the points with significant decay mode BR(t — > b)(i) with very small m-+ — m^o, 
as we have mentioned before. 

It should be noted that the BR(t — > bxf) can not simply be determined by m^f- When 
the mass of LSP m^o is much smaller than that of lighter chargino ^^+, the decay mode 
t — > txl is still significant. Furthermore, if the kinematics is allowed, other decay modes 
t — > txl an d t — > bxt are open and can reduce the decay mode BR(t — >• bxl)- For both 
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cases, the constraints from this search become weaker or even invalid. 
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TABLE IV: Summary of cuts and observed 95% CL upper limits on the excess event number, 
following the ATLAS b-jets+MET analysis for 2.05 fb" 1 Q- 
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FIG. 14: The ratio of tt*+bb* event number and observed 95% CL upper limit versus stop mass in 
the signal regions SR0-B2 and SR0-C2. The color scale indicates neutralino mass. 



Next, we consider the ATLAS b-jets+MET search based on 2.05 fb" 1 of data 751 ] . In this 
analysis, the number of b-jets is required to be at least one or two. Moreover, two signal 
regions allow one lepton in the final states. This search can be supposed to constrain the 
the signal pp — > t\t\. The corresponding cut conditions and limits are summarized in Table. 
fllVI) . The ratios N/Nn m from two most stringent channels SR0B2 and SR0C2 are shown in 
Fig. (jT4~I) . We can see even these channels require a large m e ff cut as m e fj > 700 — 900 
GeV, they can still exclude many parameter points. It is remarkable that the channel 
SR0B2 (SR0C2) can exclude signals with maximum m~ tl up to 380 GeV(480 GeV). The 
constraints from lb-jet+MET signal regions are weaker than 2b-jets+MET due to the large 
backgrounds. The lepton + b-jets + MET in the same analysis also can not achieve better 
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constraints and are omitted here. 
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TABLE V: Summary of cut conditions, following the CMS b- jets + MET analysis for 4.98 fb -1 . 
The capital letter "L" and "T" means "loose" and "tight" respectively. 
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FIG. 15: The ratio of tt*+bb* event number and observed 95% CL upper limit versus stop mass in 
the signal regions 1BL and 2BL. The color scale indicates neutralino mass. 



Then we consider CMS b-jets+MET search based on 4.98 fb -1 of data 76||. The cor- 
responding cut conditions are summarized in Table. (jV|). In this analysis, H-r defined as 
scalar sum of momenta of all the energetic jets are used to set cut conditions. A(f> norm is 
normalized azimuthal separation between f, direction and jets. The collaboration has not 
yet provided limits on the number of new physics explicitly. Here we use the formula given 
in Ref. 



to roughly estimate the upper-limits, 



X 



obs 



N SM - N, 



BSM 



^ BSM "I" N S M T u SM 

where iV & s is the number of observed events, Nsm and usm are the predicted background 
number and the uncertainty due to statistic and systematic reasons respectively. Requiring 



(13) 
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X 2 < 3.84, we can get the upper-limits on four channels as 124, 14.7, 58 and 41 respectively. 
Because this analysis dose not require very hard leading jet, the numbers of events passing 
all the cuts conditions are larger than last analysis. The ratios N/Nu m are given in Fig. [15j 
Because the "tight" cut conditions require very large $ the signals are significant reduced 
in these signal regions. Here we only show the constraints from "loose" searches. 

The bounds obtained from inclusive and exclusive 2-bjets+MET searches can be com- 
pared from Fig. ffl~2]) . Fig. f|T6|) and Fig. (JI5]) . Roughly speaking, the bounds from inclusive 
2-bjets+MET searches seem to be more stringent, as demonstrated by the fact the selected 
points are squeezed in a smaller range of N/Nu m between 0.1 and 1 and the maximum 
excluded can reach to ~ 400 GeV. 
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TABLE VI: Summary of cuts and observed 95% CL upper limits on the excess event number, 
following the ATLAS "heavy top" llepton+b-jets+MET analysis for 4.7 fb" 1 |76|. 
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FIG. 16: The ratio of tt*+bb* event number and observed 95% CL upper limit versus stop mass in 
the signal regions B and D. The color scale indicates neutralino mass. 



Since the semileptonic mode of pp — > t\t\ — > ttXiXi can have a large branching fraction 
and enjoys a smaller SM background, so it is expected the search for this mode should 
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be stringent. Then we consider the constraints from ATLAS llepton+b-jets+MET searches 



based on 4.7 fb _1 of data 77J] 5 - The corresponding cut conditions and limits are summarized 
in Table. (IVII) . It requires that the number of isolated lepton is exactly one. Obviously, this 
analysis is optimized for searching stop pair production with decay mode t — > t\\- One top 
from stop decay is required to decay hadronically and the other semileptonically. rrix is the 
transverse mass defined as 



m T = ^2p l T ft(l - cosA0(7,|y). (14) 

rriT denotes the mass scale of mother particles which decay into charged leptons, and my 
cut can be used to reduce W + jets backgrounds. To suppress the backgrounds from dilep- 
tonically decaying top pair, a specific cut on three-jet invariant mass rrijjj is required. Two 
jets with rrijj > 60 GeV and smallest AR are assumed to be originated from a hadronically 
decaying W boson, and a third jet which is closest to the reconstructed W boson is selected. 
These three jets may be the decay products of a hadronically decaying top, and the invariant 
mass rrijjj is required to be around top mass 130 GeV< rrijjj < 205 GeV. 

The ratio N/Nn m is shown in Fig. (TT6]) . It is obvious that this channel is sensitive to the 
mass range of stop from 270 GeV to 400 GeV, when the main decay mode of stop can be 
t txi- On the other hand, for the parameter points with heavy stop > 400 GeV and light 
LSP < 150 GeV in our scan, the main decay mode of stop t — > tx®, t — > fx® or t — > bxt are 
open. In this case, although the mass difference between stop and LSP is large, the limits 
becomes weaker due to the absence of energetic top quarks in the final states. From Fig. 
(]16p we can also see the constraints from signal region B are more strong than region E due 
to the smaller ft cut condition. 

In Fig. (TTTj) . we summarize the most stringent bounds from direct SUSY searches for 
both stop and sbottom pair productions. It is observed that for the stop pair production, 
the current LHC collaborations can exclude stop(sbottom) up to 400 GeV or so when only 
the exclusive signals are considered. While when the inclusive signals are considered, the 
direct searches can exclude the mass of stop up to 500 GeV, as shown by the bottom right 
plot. 



5 The ATLAS analysis of lepton+b-jets+MET in Ref. 78[ is also optimized for stop pair production with 
decay mode t — > txi, while tops are assumed to decay hadronically. The constraints from this channel are 
weaker than the semileptonic channel over most of the parameter space. 
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FIG. 17: Top left: the maximum value of N(tt*)/Nn m (which is determined by considering all 
N(tt*)/Nn m from different channels) versus stop mass, where the color scale indicates chargino 
mass. Top right: the maximum value of N(bb*)/Nn m versus sbottom mass, where the color scale 
indicates neutralino. Bottom left: the maximum value of (N (tt*) + N (bb*)) /Nn m versus stop mass, 
where the color scale indicates neutralino mass. Bottom right: the maximum value of N(bb*)/Nn m 
versus the maximum value of N(tt*)/Nu m , where the color scale indicates stop mass. 



Finally, it is should be mentioned that CMS and ATLAS collaboration have released a 

:h 5 fb" 1 of data 



series of results of light stop/sbottom searches wit 
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80h83| . Some results 



84|. In these searches, some powerful 
and "Razor" 81] . are 



based on y/s = 8 TeV have also been provided 
methods based on additional kinematic variables, such as " ar ' 
performed. These analyses can set very stringent constraints on the sbottom and stop 
pair production with large BR(b — > an d BR(t —> tx?)- Some parameter points with 
mi > 500 GeV may also be excluded. Here we do not perform these analyses, and leave 
them in the future studies. 
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C. benchmark points 



Point 
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A K (GeV) 
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-785 


Mq 3 (GeV) 
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M in = M ln ( GeV ) 


998 
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M h = M i R (GeV) 
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140 


Mi (GeV) 
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M 2 (GeV) 
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160 



TABLE VII: Input parameters of four bench mark points. 



In this subsection we choose four benchmark points that have passed all the constraints 
we considered above and discuss their features at the LHC. In Table . fIVIII) and Fig. f|T8|) . 
we tabulate the mass spectra as well as main decay modes of stop and sbottom of four 
benchmark points. In the first and fourth bench mark point, the SM-like Higgs boson is 
Hi\ while in the second and third bench mark point, the SM-like Higgs boson is H\. For 
all four bench mark points, the Hd dominated Higgs bosons, including H%, A2 and H ± are 
quite heavy and are decoupled. 

It is also noticed that the lightest stop and sbottom are left handed dominant in the 
first two bench mark points while are right-handed dominant in the other two bench mark 
points. Due to the branching fractions BR(ti — > t\\) of these bench mark points are less 
than 50% and also partially due to the mass splitting between bm~ t = — m^o , such four 
bench mark points have not been excluded by LHC searches. 

There are some interesting phenomenologies for these benchmark points. It is noticed 
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20.7% 


-> txr) 


82.9% 
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TABLE VIII: Mass spectra as well as main decay modes of stop and sbottom in four bench mark 
points, where the label W* in the row marked by Br(xf —> W^Xi) means the off-shell W boson. 
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that the branching fraction of ti — > tx 2 an d t\ ^xH can be quite large for the first two 
benchmark points. For the first benchmark point, the x 2 can dominantly decay to H 2 x\ with 
branching fraction 76.4%, then the signature pp — > t\t\ — > tiH 2 H 2 XiXi can be sizable when 
luminosity is large enough. Therefore multi-bjets plus top pair searches should be useful for 
this channel. For the second benchmark point, the x\ can dominantly decay to x? phis an 
off-shell Z boson with branching fraction 100%, then the signature pp — > t\t\ — > tiZ*Z*XiXi 
(W*, Z* means off-shell gauge bosons) can be sizable when luminosity is large enough. 

For the first benchmark point, the X3 can dominantly decay to Z(Ai)xi with branching 
fraction 66.2% (31.3%), then the signature pp — > tit* — > tiZ(Ai)Z(Ai)xiXi can be sizable. 
For the second benchmark point, the X3 a lso goes to Z(Ai)xi but with branching fraction 
83.0% (16.9%). Therefore, except the pp —> tixxXi search, the search for pp — > t\i\ — >• 
tiZ(Ai)Z(Ai)x < iXi can be complementary to constrain these two benchmark points (here 
the branching fraction of A\ — > bb can be 90% for both benchmark points). 

It is also remarkable that due to the fact that the lighter sbottom quarks are dominantly 
left-handed, consequently its dominant decay mode is b\ — > tx{ with branching fraction of 
82.9% and 93.2% respectively. For both benchmark points the chargino xt decays 100% to 
Xi and W boson. Therefore the production of pp — > t\t\ can leads to a sizable final state 

with bbW + W-W + W-XiX°i- 

Obviously, the same sign lepton plus jets signature can help to constrain these two bench- 
mark points. With more dataset accumulated at \/s = 8 TeV [3], it is expected that the 
direct searches of LHC can either rule out or discover these three benchmark points. 

For the third benchmark point, the decay mode t± — > bxt is quite large, while xt dom- 
inantly goes to Zxt and VF + Xi with branching fraction 44.4% and 41.3%, respectively. 
Therefore the production of pp — > tit* can leads to sizable final states like bbW + W~x\x\- 
The decay mode b\ — > bx% is the dominated b\ decay mode, while X2 dominantly goes to 
W*xi- 

For the forth benchmark points, more decay modes are open and many decay modes share 
large branching fractions. In the ti — > bxt m ode, the heavier chargino xt dominantly goes 
to H 2 /Hi/Zxf with branching fractions 36.0%/10.8%/14.7% and Wx\/x\ with branching 
fractions 15.0%/10.5%. For the sbottom, the decay chains also become quite long. It might 
be more challenging to exclude or discover this benchmark point. 

For comparison, we tested the benchmark point 3 given in [3] and found it is still alive. 
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FIG. 18: Mass spectra of benchmark point 1 (top left), benchmark point 2 (top right), benchmark 
point 3 (bottom left ), and benchmark point 4 (bottom right). 

The main decay modes of light stop in this point are t\ — >■ tx\, t\ — > and t\ — > bxt 
with branching fractions 16%, 29% and 45%, while x[> an d xt decay to Z*Xi and Z/H 2 xt 
with large branching fractions respectively. We also tested the benchmark points given in 
[l^ and found they can survived from the search of the direct LHC searches. The main 
decay modes of light stop in these two benchmark points are ti — > t\\ and t\ — > bxf with 
branching fractions 23% and 56% (20% and 68%). In these benchmark points, compared 
with the MSSM case, the branching fractions of b\ — > bxi and t\ — > tx\ can be suppressed 
drastically if Xi is singlino dominant. Such a fact leads to weaker constraints when we apply 
the LHC direst searches to parameter spaces with lighter stop and sbottom. 
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IV. CONCLUSIONS 



The NMSSM provides a natural framework for the recently discovered 125 GeV Higgs 
boson. Within the NMSSM, we have analyzed the constraints from the 125 GeV Higgs boson 
as well as the results from the dark matter searches to the parameter space. We concentrate 
on the LHC direct SUSY searches on the allowed parameter points. 

We have focused on scenario where the stop/sbottom can be lighter than 500 GeV and 
performed a detailed study to examine how the SUSY direct search can constrain them by 
using the results based on y/s =7 TeV and 2 ~ 5 fb _1 of data. It is found that the direct 
SUSY searches, especially the channels with tagged b jets, are powerful and can put bounds 
to the allowed parameter space. 

We would like to point out that when the inclusive signatures of both stop and sbottom 
pair productions are considered, the direct SUSY searches can exclude many parameter 
points with the left-handed stop/sbottom up to 500 GeV or so. With y/s =8 TeV and 5 
fb _1 of data or more, although kinematics could be a little different, we can expect the direct 
SUSY searches will push the light third generation squarks of NMSSM to narrower corners. 
However, for benchmark points given in our work and Ref. 12J, LL3J], special strategies and 
kinematic variables are still needed for searching light stop/sbottom pair signatures. 
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